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ABSTRACT 

We present a catalog of 191 extragalactic sources detected by the Atacama Cosmology Tele- 
scope (ACT) at 148 GHz and/or 218 GHz in the 2008 Southern survey. Flux densities span 14 
to 1700mJy, and we use source spectral indices derived using ACT-only data to divide our 
sources into two subpopulations: 167 radio galaxies powered by central active galactic nuclei 
(AGN), and 24 dusty star-forming galaxies (DSFGs). We cross-identify 97% of our sources 
(166 of the AGN and 19 of the DSFGs) with those in currently available catalogs. When 
combined with flux densities from the Australian Telescope 20 GHz survey and follow-up 
observations with the Australia Telescope Compact Array, the synchrotron-dominated popu- 
lation is seen to exhibit a steepening of the slope of the spectral energy distribution from 20 
to 148 GHz, with the trend continuing to 218 GHz. The ACT dust-dominated source popula- 
tion has a median spectral index, 0 : 148-2185 °f 3.7lg'J 6 , an ^ includes both local galaxies and 
sources with redshift around 6. Dusty sources with no counterpart in existing catalogs likely 
belong to a recently discovered subpopulation of DSFGs lensed by foreground galaxies or 
galaxy groups. 

Key words: galaxies: surveys - galaxies: active - millimeter: galaxies 
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1 INTRODUCTION 

The technologies that enable observations of large numbers of 
millimeter and submillimeter sources were developed relatively re- 
cently. They open up for study a previously unexplored regime that 
has the power to reveal the evolution of underlying galaxy popula- 
tions over cosmic time, and in particular over epochs of intense 
star formation. Instruments such as the Submillimeter Common 
User Bolometer Array (SCUBA, SCUBA-2; Holland et al. 1999, 
2013) operating at 850 /mi, the Balloon-borne Large Area Submil- 
limeter Telescope (BLAST; Pascale et al. 2008) operating at 250, 
350 and 500 /tm, the Large Apex Bolometer Camera (LABOCA; 
Siringo et al. 2009) operating at 870 /tm, and the AzTEC millime- 
ter wavelength camera (Wilson et al. 2008) operating at 1.1 and 
2.1 mm have mapped up to 10 deg 2 of the sky, but greater coverage 
has been limited by the large amount of integration time required 
to conduct blind surveys to significant cosmological depth. In- 
creased sky coverage and sensitivity motivated the construction of 
space-based observatories such as Spitzer (Werner et al. 2004) and 
Herschel (Pilbratt et al. 2010), operating in the wavelength regime 
3-500 /an and covering up to tens of square degrees. At longer 
wavelengths, the Wilkinson Microwave Anisotropy Probe ( WMAP\ 
Wright et al. 2009) covered the whole sky at frequencies up to 
94 GHz (A = 3.2 mm), but observations of unresolved extragalac- 
tic sources were only complete above 2Jy due to its large beam 
size (roughly 13' at 94 GHz). The Planck space telescope con- 
tains a high frequency instrument observing at frequencies span- 
ning 100-857 GHz (Lamarre et al. 2010). At 143 and 217 GHz (2.1 
and 1 .4 mm), the Planck beam sizes are approximately l' and 5' re- 
spectively, allowing for a source catalog that is complete down to 
flux densities of about 1 Jy. Large-area (> 100 deg 2 ) ground-based 
radio surveys probe with higher resolution than WMAP or Planck, 
but only as high in frequency as 20 GHz (e.g., the Australia Tele- 
scope 20 GHz survey; Murphy et al. 2010). Thus there is a niche 
for millimeter wavelength large-scale mapping experiments with 
spatial resolution superior to that of the space-based observatories. 

One such experiment is the Atacama Cosmology Telescope 
(ACT; Swetzetal. 2011). The ACT collaboration released Cos- 
mic Microwave Backgound (CMB) temperature anisotropy maps 
made with arcminute resolution from its 2008 observing season 
at 148 and 218 GHz (Diinneretal. 2013). These maps also con- 
tain galaxies that are luminous at millimeter wavelengths. Mea- 
surements of the millimeter fluxes of a large sample of sources 
have great potential to discriminate among source population 
models and reveal new source populations. For example, recent 
millimeter surveys (e.g., Vieira et al. 2010; Marriage et al. 2011; 
Planck Collaboration VII 2013) have yielded source counts that led 
to updated source models (Tucci et al. 2011) and the discovery of a 
new class of high-flux dusty galaxies (this publication; Vieira et al. 
2010; Negrello et al. 2010). 

ACT-detected sources at flux densities greater than 20 mJy are 
predominantly blazars that have synchrotron-dominated spectral 
energy distributions (SEDs). These galaxies are powered by active 
galactic nuclei (AGN), central black holes that accrete nearby mat- 
ter in a process that produces time-variable jets out of the plane of 
the accretion disk; in the case of blazars, the jet direction lies close 
to our line of sight. Ejected ionized particles spiraling around mag- 
netic field lines in the jets create the observed synchrotron emis- 
sion. The study of blazar SEDs and source counts provides infor- 
mation about AGN physics (e.g., de Zotti et al. 2010; Tucci et al. 
2011; Planck Collaboration XV 2011). 

The second population identified at ACT wavelengths is com- 


prised of infrared-luminous, dusty star-forming galaxies (DSFGs), 
which exhibit modified blackbody emission at sub-millimeter to 
millimeter wavelengths, diminishing toward longer wavelengths. 
The observed SED is dominated by thermal emission from dust 
grains that have been heated by the prodigious optical and ultravi- 
olet flux produced by newly formed stars (e.g., Draine 2003, and 
references therein). A fraction of these galaxies may have a signifi- 
cant non-thermal contribution from AGN at their cores, in addition 
to the thermal dust emission. 

Some of these DSFGs are local galaxies (z<Cl), cross- 
identified with Infrared Astronomy Satellite (IRAS) sources 
(Devereux & Young 1990). The IRAS wavebands (A= 12-100 pm) 
are relatively insensitive to dust with temperatures below 30 K, 
a significant and largely unexplored component of many nearby 
galaxies (Planck Collaboration XVI 2011). Recent results from 
Herschel (Amblard et al. 2010) and BLAST (Wiebe et al. 2009) 
have begun to extend our picture of the cold dust in galaxies, but 
millimeter-wavelength experiments can, through probing dust in 
nearby galaxies, contribute to establishing a well-calibrated dust 
SED for typical, low-redshift DSFGs. This in turn has cosmologi- 
cal ramifications, as current analyses of the first generations of stars 
and galaxies that fuel the cosmic infrared background (CIB) rely on 
understanding and extrapolating from dust SED templates. 

High-redshift DSFGs were observed by SCUBA in the first 
systematic survey of these sources, which create a significant 
fraction of the CIB emission (Blain et al. 2002). Selected at 
850 /tm, these sources came to be known as submillimeter galax- 
ies (SMGs). This new population of galaxies subsequently became 
a focus of observations (e.g., WeiB et al. 2009; Viero et al. 2009; 
Austermann et al. 2010) and stacking analyses that resolved more 
of the CIB into emission from discrete, dusty star-forming galax- 
ies (e.g., Dole et al. 2006; Devlin et al. 2009). Though most SMGs 
will be undetectable by current millimeter survey instruments, a 
new population of sources that are significantly brighter and rarer 
than the submillimeter-selected SMGs and that similarly exhibit 
dust-dominated spectral indices has been identified at millimeter 
wavelengths (Vieira et al. 2010). They do not have counterparts in 
the IRAS catalog, indicating that they are not members of the stan- 
dard, local, ultraluminous infrared galaxy (ULIRG) population. 

These sources have recently been shown to be a new, 
higher redshift (z > 3) subpopulation of the progenitor galaxy 
background, brought to the fore in millimeter wavelength sur- 
veys because they are lensed by foreground galaxies or galaxy 
groups (Negrello et al. 2007; Lima et al. 2010; Negrello et al. 
2010; Vieira et al. 2013; Hezavehetal. 2013). Such objects are 
extremely rare (at submillimeter wavelengths, for example, see 
Rex et al. 2010; Lupu et al. 2012), but wide area surveys will find 
more. The South Pole Telescope (SPT; Carlstromet al. 2011) has 
reported significant numbers of these sources. This population pro- 
vides an avenue for follow-up research to study the details not only 
of lensed SMGs, but also of the lens systems (e.g., Ikarashi et al. 
2011; Scott et al. 2011; Lupu et al. 2012; WeiB et al. 2013). 

In this paper, we report for the first time the discovery of 
DSFGs in the ACT data. This is the first multifrequency analysis 
of the ACT sources and the second report on ACT extragalactic 
sources. Marriage et al. (201 1), hereafter Mil, presented a catalog 
of sources at 148 GHz only from the 455 deg 2 of the 2008 ACT 
Southern survey with the best uniformity and coverage. Here we 
extend those results to include 218 GHz (and updated 148 GHz) 
flux densities for the same 2008 ACT Southern survey region from 
the new (Diinneretal. 2013, hereafter D13) data release maps. 
Gralla et al., in prep, will present sources detected in the ACT 
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Figure 1 . Sensitivity maps with source detections at 148 GHz (top) and 218 GHz (bottom), showing the most uniform 455 deg 2 patch of the Southern strip with 
the greatest depth of coverage. It lies between right ascension 00 1 1 2" 1 and 07 h 07 m , and declination —56° 11' and — 49°00'. The deepest data correspond to 
an exposure time of 23.5 minutes per square arcminute and a lcr sensitivity of 2.34 mJy at 148 GHz and 3.66 mjy at 218 GHz. White circles or squares mark 
the locations of ACT sources with a size proportional to the log of the associated source flux density. Circles denote sources designated as AGN, and squares 
denote sources with spectra indicative of DSFGs. Toward the edges of the map, the variation in local noise properties due to uneven coverage is more apparent. 


Equatorial strip. Mapmaking of the ACT 277 GHz dataset is cur- 
rently under way. 

The layout of the paper is as follows. Section 2 describes the 
2008 season ACT observations and the reduction of raw data into 
maps, as well as follow-up observations made with the Australian 
Telescope Compact Array (ATCA). Section 3 details our method of 
source extraction. The source catalog, including its astrometric and 
flux density accuracy, and its estimated completeness and purity, is 
discussed in Section 4. Section 5 compares our catalog with cur- 
rently available datasets; Section 6 gives the source number counts. 
Trends observed in the spectral indices of our source populations 
are analyzed in Section 7. We conclude in Section 8 with a sum- 
mary of our results. 


2 OBSERVATIONS AND DATA 
2.1 ACT Observations 

The ACT experiment (Swetzetal. 2011) is situated on the 
slopes of Cerro Toco in the Atacama Desert of Chile 1 at an ele- 
vation of 5,190 m. ACT’s latitude gives access to both the north- 
ern and southern celestial hemispheres. Observations occurred si- 
multaneously in three frequency bands, at 148 GHz (2.0 mm), 
218 GHz (1.4 mm) and 277 GHz (1.1 mm) with angular resolutions 
of roughly 1.4' , 1.0' and 0.9' respectively. Observations of Saturn 
were used to determine beam profiles and pointing (Hincks et al. 
2010; Hasselfield et al. 2013). From 2007 to 2010, ACT targeted 
two survey regions: the Southern strip centered around 8= —52.5° 
and the Equatorial strip centered around 5 = 0°. Further informa- 
tion about the ACT observations can be found in D13. 


1 22.9586° south latitude, 67.7875° west longitude. 


2.2 ACT Data 

The reduction of raw ACT data into maps is detailed in D13; 
we will briefly review that process here. Each ACT detector ar- 
ray (one per frequency band) is composed of 1,024 detectors, with 
each detector timestream first analyzed and then kept or rejected 
based on multiple criteria, such as telescope operation, weather 
conditions, cosmic ray hits, or other interference. Approximately 
800 (700) hours of data from the 2008 Southern strip for 148 GHz 
(218 GHz) remain after these cuts. 

Maximum likelihood maps of pixels 30" on a side are pro- 
duced from the timestream data. From an initial estimate of the 
maps, source profiles are fit for. We then subtract the timestream 
models for the point sources from the data and re-map. This pre- 
vents point source power from being aliased into the map, and im- 
proves the final flux density estimates for S/N > 5 sources (quan- 
tified below). Similarly, an estimate of the CMB signal is also re- 
moved. The mapping equation (Tegmark 1997) is then solved itera- 
tively using the preconditioned conjugate gradient method. Finally, 
the signals that had been subtracted are added back to the map. We 
conducted simulations in which a signal corresponding to sources 
of known flux were injected into raw ACT detector time streams 
before mapping. The new mapping procedure with “source subtrac- 
tion” described briefly here and more thoroughly in D13 (Section 
1 1.4) corrects a 3-5% downward bias in the recovered flux densi- 
ties. 

The measured source flux densities have converged on a sin- 
gle value by the 25 th iteration of the maps used in our study. For 
148 GHz, the D13 release map is from iteration 1000 of the map- 
maker, and we use this map both for convenience and in order that 
this study be exactly reproducible using public data. For our study 
of source flux densities at 218 GHz, we use iteration 200, which 
is more than adequate. Between the map iterations used in this 
study, the fluxes of the sources change less than 1%. One may ask 
whether the value to which the measured source flux density has 
converged is accurate. Based on end-to-end simulations in which 
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mock sources are injected into the ACT time streams, we estimate 
this accuracy at 3%. 

The 148 GHz map is calibrated in temperature through cross 
comparison of spectra over the range 400 < i < 1000 with WMAP 
as in Hajian et al. (2011) with an uncertainty of 2%. This cal- 
ibration is transferred to the 218 GHz map (2.4% uncertainty) 
through cross-calibration with 148 GHz. Hasselfield et al. (2013), 
hereafter H13, used this calibration in measurements of the tem- 
peratures of Uranus and Saturn. The resulting planet temperature 
estimates were consistent with previous analyses (Griffin & Orton 
1993; Goldin et al. 1997). To summarize the results of H13, both 
previous analyses adopted a flux standard based on a model of 
Mars (Wright 1976; Ulich 1981) with 5% systematic errors. The 
HI 3 estimate for the temperature of Uranus and Saturn was 5% be- 
low those of the previous analyses, a deficit which was also seen in 
the WMAP analysis of Weiland et al. (2011) at 94 GHz and is most 
likely due to systematic errors in the Mars standard. The agreement 
between the H13 analysis, calibrated on extended CMB emission, 
and the previous analyses, calibrated on Mars, lend confidence that 
the ACT WMAP -based calibration at large angular scales will also 
apply to the calibration of sources at small angular scales. 

Photometry based on matched filtering (Section 3) relies on 
an accurate estimate of the instrument beam shape. The impact 
of beam error on photometry depends on the details of the filter 
(Equation 1 ), which in turn depend on both the beam window func- 
tion and the angular power spectrum of the background and noise. 
We have found that the error in the solid angle of the beam is a 
conservative proxy for the photometric error due to beam uncer- 
tainty, and we use this simpler quantity (solid angle uncertainty) in 
the beam-related photometry error estimates described here. Due 
to uncertainties in the beam measurements from Saturn we as- 
sign 1% photometric errors for both bands. The profiles of bright 
point sources in the survey maps suggest that the effective beam 
for the survey is slightly broader than that measured from Saturn. 
The survey maps consist of overlapping observations taken over 
the entire season, and the observed broadening may be attributed 
to a night-by-night jitter in the telescope pointing with an rms of 
5 ± l"(Hasselfield et al. 2013). The uncertainty in jitter correction 
corresponds to 1% (148 GHz) and 2% (218 GHz) photometric un- 
certainties, which are correlated between bands. 

Finally, for our photometric uncertainty budget we account for 
the fact that AGN have a lower effective frequency band center 
(and thus slightly broader beam) than Saturn, and DSFGs have a 
higher effective frequency band center (and thus narrower beam; 
Swetzetal. 2011, Table 4). We choose effective 148 GHz and 
218 GHz frequency centers corresponding to halfway between a 
steep spectrum AGN and a DSFG: 148.65 GHz and 218.6 GHz. 
This choice introduces a photometric bias of less than 1.5% at 
148 GHz and 1.1% at 218GHz, which is positive for DSFGs and 
negative for steep spectrum AGN. We fold this photometric bias 
from the source spectrum in with uncertainties due to mapping 
(3%), WMAP calibration (2%, 2.4%) and beam shape (1.4%, 2.2%) 
to obtain an overall flux density calibration uncertainty of 4.1% at 
148 GHz and 4.6% at 218 GHz. As shown in Section 5, the ACT 
flux densities agree with independent measurements to within this 
margin of error. 

For this study, we have used the most uniform 455 deg 2 of the 
2008 ACT Southern strip at 148 and 218 GHz. This region, shown 
in Figure 1, spans declination -56.2° < <5 < -49.0° and right ascen- 
sion 00 h 12 m < a < 07 h 07 m . At 148 GHz we use the data publicly 


released with D13 2 . Typical white noise levels in this region of 
the map are 30 /rK-arcminute at 148 GHz and 50 /rK-arcminute for 
218 GHz. As described in Section 3, when matched-filtered with 
the ACT beam, this white noise level results in a 1 a point source 
flux density sensitivity in the best covered regions of <to = 2.34 mJy 
at 148 GHz and 3.66 mJy at 218 GHz. The sensitivity levels in Fig- 
ure 1 are proportional to the square root of the number of obser- 
vations at that map location, N 0 b s , with one observation per 0.005 
seconds. Then the sensitivity level in a given portion of the map is 
rr oo \J Wobs ,max /Nobs- 


2.3 ATCA Observations 

Marriage et al. (2011) found that the sample of ACT 148 GHz- 
detected sources cross-identified with the Australian Telescope 
20 GHz survey (AT20G; Murphy et al. 2010) is dominated by 
sources with peaked or falling SEDs using flux densities mea- 
sured at 5, 20 and 148 GHz. The study also confirmed the find- 
ings of the AT20G study (Murphy et al. 2010; Massardi, M. et al. 
2011), namely that this population of radio sources is character- 
ized, on average, by spectral steepening between 20-30 GHz and 
148 GHz. However, this sample from Mil was incomplete, biased 
in a way that favored sources with negative spectral indices between 
20 and 148 GHz due to the AT20G survey completeness level of 
78% above 50 mJy at 20 GHz. ACT-selected sources with 148 GHz 
flux densities less than 50 mJy and flat or rising spectra may not 
have been detected by AT20G. Therefore, in order to complete the 
Mil 20-148 GHz spectral study, a targeted set of measurements of 
flux densities at 20 GHz was made for the Mil sources that were 
not identified with any source in the AT20G catalog within a 30" 
search radius. 

Scheduling and weather constraints permitted us to observe 41 
of these 48 sources with the 6 x 22 m antenna array of ATCA 3 over 
6.5 hours at 20 GHz on November 10, 2010, when the array was in 
its East-West 750A configuration, with 15 baselines ranging from 
77 to 3750 m. The primary beam of each array telescope at 20 GHz 
is 2.3' , with a resolution of 0.5" possible with the full array in this 
configuration. The synchrotron-dominated spectra of these sources, 
as revealed in Mil, indicated that they were compact (AGN), and 
would be unresolved. For these observations, the two 2 GHz-wide 
frequency bands of the new ATCA Compact Array Broadband 
Backend (CABB; Wilson et al. 2011) digital array correlators were 
set to be adjacent by centering them at 19 GHz and 21 GHz. The re- 
duced average flux densities over the whole bandwidth of the corre- 
lator corresponds to the “20 GHz” flux density. Good weather con- 
ditions prevailed throughout the observations, from 07:35:24.9 UT 
through 14:1 1:54.9 UT. A few data blocks were flagged and re- 
moved in order to minimize noise and any spurious effects. 

The primary flux calibrator used for these observations was 
PKS B 1934-638, and the bandpass calibrator used was PKS 
B1921-293. The target sources were expected to have flux densi- 
ties at or just below the AT20G survey limit of 40 mJy. Each target 
source was observed once for 1.5 minutes for an rms noise level of 
^ 0.15 mJy/beam. Target observations were interleaved with obser- 
vations of four secondary calibrators for pointing and phase correc- 
tions, chosen to lie close to the targets. 


2 http://lambda.gsfc.nasa.gov/product/act/ 

3 http://www.narrabri.atnf.csiro.au/ 
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2.4 ATCA Data 

The ATCA follow-up data was reduced using a fully au- 
tomated, custom, shell script pipeline based on tasks from the 
MIRIAD aperture synthesis reduction package (Saultetal. 1995). 
An initial inspection of the data was performed to identify contam- 
ination or any problems in the data acquisition. Automatic proce- 
dures were used to identify and flag data for each frequency band 
affected by shadowing or known radio contamination, resulting in 
less than 1% of the band being flagged and cut. The pipeline then 
generated the calibration solutions for bandpass, flux density ampli- 
tude and phase based on the calibrators. After checking that results 
were consistent between the two frequency bands, we merged the 
calibrated target visibilities before extracting the flux densities to 
improve sensitivity. 

Analyses of the AT20G survey (Murphy et al. 2010) and of 
other ATCA projects (Massardi etal. 2011; Bonavera et al. 2011) 
have shown that the triple correlation technique (Thompson et al. 
1986) can be effectively used to determine source flux densities 
down to tenth of mjy scales for point-like sources unaffected by 
poor phase stability. This technique is effective for point sources 
even if they are not in the phase center and in the case of low S/N. 

We manually inspected the 20 GHz deconvolved images, and 
wherever a source was clearly identifiable we measured its flux 
density and flux density error at its peak in the image and the rms 
of the image pixels in a region unaffected by the source emission. 
For sources with high S/N, deconvolution techniques were able to 
reconstruct images, despite the poor uu-coverage of our observa- 
tions. In the few cases where the source appeared extended with 
respect to the synthesized beam we estimated the flux density by 
integrating over the source area. In some cases, the image showed 
a confused field with multiple peaks, primarily due to sidelobes; 
for these cases we extracted a flux density estimate using the triple 
correlation technique. The flux density errors were obtained from 
the rms of the visibilities of the V Stokes parameter. This assumes 
that these objects have negligible circular polarization, which is an 
overestimate in the cases where the source is not in the phase cen- 
ter. The same statistic (V Stokes mis) was used to indicate the noise 
level reached for fields where no source was identified at 20 GHz. 
Of the 41 ACT sources observed in the follow-up campaign de- 
scribed here, all had a measured flux density at 20 GHz. The final 
results of this study are included as part of Table 4 with asterisks 
indicating sources that were part of the follow-up ATCA sample. 


3 ACT SOURCE EXTRACTION 

We use the same matched filtering method as described in Ml 1 
to produce an estimate of the amplitude of a source at the location 
of the source center. In this section we summarize the method with 
an emphasis on aspects for which this analysis differs from Mil. 

Initially, the data are weighted by the square root of the num- 
ber of observations per pixel. This method results in a map with an 
approximately constant white noise level and a natural apodization 
as the number of observations falls off toward the edge of the survey 
region. Using a new procedure, we filter the entire D13 148 GHz re- 
lease map and corresponding 218 GHz map, both of which extend 
beyond the deep 455 deg 2 area of this study, in order to eliminate 
any potential edge effects. We then use only the 455 deg 2 region 
for source identification. 

The following filter is applied to the Fourier transfomi of the 

map: 



Figure 2. Matched filters (solid lines) are band pass filters that smooth the 
map on the scale of the beam (dashed lines) and remove large scale structure 
associated with the CMB, other astrophysical signals (e.g., the Sunyaev- 
Zel’dovich effect) and residual contamination from atmospheric brightness 
fluctuations. The beam and filter functions are plotted with unit normaliza- 
tion at their peak (9 = 0). The matched filters were binned in radius to better 
show the relevant angular scales. In the analysis, we use two dimensional 
filters to capture the anisotropic character of the noise. 

$(k) = Fk 0 ,k x (k)H*(k)| T other (k) |~ 2 

f B*(k')Fk 0 ,k x (k')| T 0 ther(k') | _ 2 B(k') dk' ’ 

where k = (k x , k y ) is the angular wave number, and x and y re- 
fer to the right ascension and declination directions. B(k) is the 
Fourier transfomi of the “effective” instrument beam. As described 
in Section 2, the “instantaneous” beam is derived from planet obser- 
vations (Hasselfield et al. 2013). The “effective beam” in the 2008 
survey map is broadened by imperfect telescope repointing with 
mean deviation ag = 5". This broadening is included in B( k) by 
multiplying the instantaneous beam transform (released with D13 
for 148 GHz) by exp(— £ 2 ag/2). T’other is the Fourier transform 
of all components of the data besides point sources (i.e., atmo- 
spheric or detector noise, CMB, etc.). The function Fk 0: k x ( k) is 
a high-pass filter that removes undersampled large scale modes be- 
low ko = 1000 and modes with | k x \ < 100, which are contaminated 
in a fraction of ACT data by telescope scan-synchronous noise. 
While -B(k) is well approximated as azimuthally symmetric, we 
retain the full two dimensional power spectrum ITotherjk') | 2 to 
downweight anisotropic noise in the maps. The azimuthally-binned 
1 48 and 218 GHz real space filters and associated beam profiles are 
shown in Figure 2. 

The power spectrum | Tothe^(k , ) | 2 used in this analysis was 
constructed in a different manner from Mil: we used the power 
spectrum of the data itself instead of the average of difference 
(noise) maps and models for the CMB and other contaminating sky 
emission. This estimate is robust since the total power from the ex- 
tragalactic source signal is low compared to the CMB, atmospheric 
noise, and white noise. To avoid a noisy estimate of the power spec- 
trum, we smooth the power spectrum | T ot her(k') | 2 with a Gaus- 
sian. The exact fonnulation of the smoothing does not significantly 
change the resulting filtered map. 

Applying the matched filter can cause ringing in the maps 
around the very brightest sources, which impacts source extraction 
around other, low S/N, sources. Therefore, we identify sources with 
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Figure 3. Filtered 148 GHz (top) and 218 GHz (bottom) submaps. The data 
have been match-filtered such that the grey-scale is in units of flux density 
(mJy) with white (black) corresponding to -lOmJy (30mJy). Insets show 
the flux density distribution across the data as a grey histogram that has a 
standard deviation of2.57mJy (3.78 mJy) for 148 GHz (218GHz). Several 
sources, both synchrotron and dust dominated, marked as circle and square 
outlines respectively, are apparent as black beam-sized flux excesses. The 
white extended object in the 148 GHz map, marked by a diamond, is the 
Sunyaev ZeTdovich effect decrement from Abell 3128 NE (ACT-CL J0330- 
5228). 


S/N > 50 in an initial application of the filter, and mask them in the 
maps. Sources with S/N < 50 are then extracted by match filter- 
ing these masked maps. Once filtered, groups of map pixels with 
S/N > 4.8 are identified as candidate sources. 

Finite map pixel size affects the measured flux densities, since 
it is rare that a source falls exactly in the center of a pixel, lead- 
ing to a systematic negative bias and increased scatter. We increase 
the map pixel resolution by a factor of 16 in a region 0.03° on a 
side around each source through Fourier space zero-padding (e.g., 
Press et al. 2007) to allow for a more precise determination of the 
source peak location, and therefore its flux density. Properly cen- 
tering the detection has the effect of boosting the S/N of typical 
S/N = 4.8 sources to ^ 5. Lastly, to account for the convolution of 
sources with the map pixel, which acts as a low pass filter (the 
square pixels become a sine function in Fourier space), we decon- 
volve the map pixel window function in the higher resolution map. 
We find that this method reduces systematic errors associated with 
pixelization to below 1%. Purity tests (Section 4.3) reveal a signifi- 
cant number of false detections below S/N of 5.0, due to local noise 
and striping artifacts. Therefore we impose a S/N = 5.0 threshold 
on sources detected in each map assuming no prior knowledge of 
source location from the other frequency. 

Once a S/N >5.0 catalog has been generated for each of the 
filtered 148 GFlz and 218 GFlz maps, we look for sources that have 
been identified at one frequency but not the other. The flux density 
for the source at the second frequency is then measured from the 
filtered map for the second frequency. We use the same prescription 
as for the flux measurement method just described, using a submap 
at the second frequency centered on the location of the detection as 
determined from the map at the first frequency. 


4 THE 148 AND 218 GHz CATALOG 

The ACT-detected source catalog is given in Table 4. We 
find 169 sources selected at 148 GHz with S/N > 5, spanning two 
decades in flux density, from 14 to 1700 mJy. The 218 GHz map in- 
dependently yielded 133 sources with S/N > 5. The combination of 
these two independent source lists gives a total count of 191, with 
1 10 galaxies detected with S/N > 5 at both frequencies. 

The catalog provides the IAU name, celestial coordinates 
(J2000), S/N and flux density estimation of each 148 and 218 GHz 
ACT-detected source. Raw flux densities are estimated directly 
from the map as described in Section 3. Deboosted flux densi- 
ties as derived according to Section 4.2 are given with associated 
68% confidence intervals. The raw and deboosted spectral indices, 
ai 48 — 2 i 8 , between 148 and 218 GHz for each source, are provided. 
The classification of sources as dust or synchrotron-dominated is 
based on the ai 48 - 2 is spectral index criterion described in Sec- 
tion 4.2. If the source was cross-identified with an AT20G catalog 
source (see Section 5), the AT20G source ID is given. If, instead, 
the 20 GHz source flux density was measured during the November 
2010 follow-up campaign (Section 2.3), the 20 GHz ID name has 
an asterisk next to it. Sources not cross-identified with one of the 
catalogs listed in Section 5 are marked with a “d” superscript. 

Correlating 148 GHz flux densities between this catalog and 
the catalog given in Mil, we find an average agreement at the 2% 
level, with larger scatter for individual sources. This level of consis- 
tency is expected given changes in calibration, map-making, beam 
profile estimates, and deboosting procedure made for this updated 
study. 

The reported flux densities are the average over approximately 
two months of observation for each source, many of which are 
AGN-driven radio galaxies and thus likely to have varied in that 
time. For example, of the three bright sources cross-identified with 
Planck, in the year between ACT and Planck observations one 
source varied by 20% in flux density whereas another did not vary 
at all within errors. However, our simultaneous multifrequency ob- 
servations allow for a consistent internal spectral characterization 
between ACT bands. In a future study, multiple years of data will 
be used to quantify the effects of variability on individual source 
spectra. 

The following sections provide details and context for the cat- 
alog values. 


4.1 Astrometric Accuracy 

Radio interferometers can achieve very precise positional ac- 
curacy for sources, so ACT-selected sources cross-identified with 
a robust radio catalog give a good measure of the positional accu- 
racy of the ACT source detections. The AT20G catalog covers the 
Southern sky, and through pointing checks against Very Long Base- 
line Interferometer (VLBI) measurements of International Celestial 
Reference Frame calibrators, the positional accuracy of AT20G is 
shown to be accurate to better than (Murphy et al. 2010). We 
exclude nearby extended/resolved sources from this analysis, deter- 
mined from cross-identification of our sources with currently avail- 
able catalogs, as distant point-like sources will present a clearer 
picture of the overall accuracy of our pointing. 

We compared the positions of the 34 ACT sources with 
S/N 148 > 16 to positions of associated sources in the AT20G cat- 
alog using a search radius of half the ACT 148 GHz beam (0.7 / ). 
Figure 4 shows the offsets in location between the AT20G right as- 
cension (RA) and declination (Dec) and the ACT-derived positions 
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Table 1. Astrometric Pointing Accuracy 



148 GHz 

218 GHz 

Mean RA Offset (") 

0.6 ±0.4 

0.5 ±0.6 

Mean Dec Offset (") 

-0.4 ±0.3 

0.0 ± 0.6 

RMS in RA Offset (") 

2.1 

3.5 

RMS in Dec Offset (") 

1.8 

3.4 



Figure 4. Astrometric accuracy of the ACT source detections. The 34 filled 
black circles are the positional offsets of ACT sources with S/N 148 >16 
from counterparts in the AT20G catalog as calculated from the 148 GHz 
map. The turquoise triangles are for the same sources, but using positions 
derived from the 2 1 8 GHz map, which has more noise. The large errorbars 
are centered on the mean offsets, and extend as far as the rms of the offsets. 
See Table 1 . 

for the sources as they were detected in the 148 GHz (black points) 
and 218 GHz (blue points) maps. The results are also summarized 
in Table 1 . The large overlaid errorbars are centered on the mean 
offsets, and extend as far as the rms of the offsets. For sources with 
lower S/N, the ACT location rms with respect to AT20G position 
becomes inflated by the effect of noise in the maps. Therefore the 
catalog uses coordinates derived from the 148 GHz map since there 
is a higher level of noise present in the 218 GHz map. 

4.2 Flux Density Accuracy 

Flux densities derived from the ACT maps have systematic 
uncertainties arising from five effects: the overall calibration un- 
certainty, the mapmaker, errors in the assumed source profile, error 
in the assumed source spectrum, and flux boosting of lower signif- 
icance candidates. The first four potential sources of flux density 
error were discussed in Section 2. Here we will discuss the last 
source of flux density error. 

The differential counts of the sources selected in our sam- 
ple fall steeply with increasing flux (Table 2, Figure 6). With no 
prior information about the source flux, the most likely scenario is 
that the measured flux is the sum of a dimmer intrinsic flux and a 
positive noise fluctuation. We use the two-band Bayesian method 
developed in Crawford et al. (2010), and report the 16, 50, and 
84 percentiles (68% enclosed, equivalent to l<r) of the posterior 


flux and spectral index distributions. For the source count priors 
in this calculation we use the sum of the models of de Zotti et al. 
(2010) for radio sources (using the Tucci et al. 2011 model re- 
sults in flux differences of <0.03cr) and Bethermin et al. (2011) 
for dust-dominated sources. Following Vieira et al. (2010), we take 
a flat prior on the spectral index between —3 and 5, consistent 
with the expected range for our populations. The two-band like- 
lihood includes negligible correlation between bands and is consis- 
tent with background astrophysical emission rather than correlated 
atmospheric emission. 

The source populations in 148 and 218 GHz naturally split 
into sources having their emission dominated by synchrotron (cen- 
tered on a = -0.6) or thermal dust (centered on a = 3.7; see Fig- 
ure 7, bottom panel). We use the threshold spectral index a = 1.66 
(Vieira et al. 2010) to divide these populations in terns of their pos- 
terior spectral index populations, with P(a> 1.66) >0.5 classi- 
fied as dusty and P(a > 1 .66) < 0.5 as synchrotron dominated. The 
classification is robust within a ± 0.5 of this threshold. While cal- 
ibration uncertainty is included in all quoted fluxes, we conserva- 
tively ignore inter-band correlation of calibration errors when flux 
deboosting to avoid making assumptions about the correlation of 
the pointing jitter, chromatic effects, and mapping/flux recovery er- 
rors (Section 2.2). This moderately boosts the error on the spectral 
index toward high S/N, but does not impact our source identifica- 
tions. 

In addition to the bias from the steepness of the population, we 
also treat the fact that the source finder locates the maximum flux 
along RA and Dec, which provides an extra two degrees of free- 
dom. This can be corrected by finding the flux after two degrees of 
freedom are subtracted from the detection significance, following 
Vanderlinde et al. (2010). Note that this departs from the treatment 
in Crawford et al. (2010) and Vieira etal. (2010). Specifically, all 
sources are corrected by a factor (S/N) / \J S/N 2 — 2. For sources 
in the range 20-25 mJy, this is a correction of 0.5% at 148 GHz 
and 1.5% at 218 GHz. We have confirmed through simulations with 
synthetic sources implanted in the ACT data that this positional de- 
boosting results in unbiased flux densities. Furthermore, we have 
compared raw flux densities from the matched filter to flux den- 
sities from the ACT data at the positions of ATCA counterparts, 
when available. The latter should not be boosted due to maximiz- 
ing the flux over position in RA and Dec. As expected, we find that 
positional deboosting accounts for the ratio between the raw flux 
densities and the flux densities derived using ATCA counterpart lo- 
cations. 

As a final consistency check, we note that for the sources ob- 
served by both ACT and Planck, flux densities are consistent at 
the « 1-2% level, suggesting that the errors assumed here are very 
conservative and free from any systematic bias (Section 5.3). 


4.3 Purity and Completeness 

The number of false detections at each frequency was esti- 
mated by running the detection algorithm on an inverted (nega- 
tive temperature) map in which we masked the sources and, in the 
case of the 148 GHz map (for which the Sunyaev-Zel’dovich effect 
was non-zero), all ACT-detected and optically confirmed clusters 
of galaxies. With this approach, no spurious detections are found 
in the 148 GHz data down to a S/N of 5. In the 218 GHz data, four 
spurious detections at S/N ^ 6 were observed with raw flux densi- 
ties in the 25-31 mJy range (Table 2). 

The purity simulations appear to be consistent with our find- 
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Table 2. Number Counts, Purity, and Completeness 11 . 


Flux Range (Jy) 


148 GHz 


218 GHz 

N sync 


N 

Purity 

Completeness 

N 

Purity 

Completeness 

0.015-0.02 

23 

100.0±0.0% 

47.1±3.2% 

8 

100.0±0.0% 

10.7±2.8% 

22 

0 

0.02-0.03 

47 

100.0±0.0% 

75.3±3.4% 

42 

92.9±4.8% 

34.2±3.5% 

46 

10 

0.03-0.05 

48 

100.0±0.0% 

96.8±1.9% 

43 

97.7±2.3% 

80.2±3.0% 

48 

10 

0.05-0.09 

20 

100.0±0.0% 

99.6±0.6% 

21 

100.0±0.0% 

98.7±0.9% 

20 

1 

0.09-0.17 

13 

100.0±0.0% 

100.0±0.0% 

12 

100.0±0.0% 

100.0±0.0% 

13 

2 

0.17-0.33 

4 

100.0±0.0% 

lOO.OiO.0% 

2 

100.0±0.0% 

lOO.OiO.0% 

4 

0 

0.33-0.65 

2 

100.0±0.0% 

100.0±0.0% 

1 

100.0±0.0% 

100.03=0.0% 

2 

0 

0.65-1.39 

1 

100.0±0.0% 

100.0±0.0% 

1 

100.0±0.0% 

100.03=0.0% 

1 

0 

1.39-2.87 

1 

100.0±0.0% 

100.0±0.0% 

1 

100.0±0.0% 

100.0±0.0% 

1 

0 


“The number of sources at each frequency with S/N = 5 over 455 deg 2 . Errors on number count are simply 


Poisson. See Figure 6 for a graph of purity/completeness-corrected differential source counts. 
^Counts of synchrotron-dominated sources are taken relative to the 148 GHz flux density. 
“Counts of dust-dominated sources are binned according to their 218 GHz flux density. 


ings from source cross-identification (Section 5). At 148 GHz, only 
one 5 a source, ACT-S J023600-530237, is not cross-identified 
with other catalogs. At 218 GHz, five 5cr sources, ACT-S J024430- 
541605, ACT-S J035034-524801, ACT-S J062747-512614, ACT- 
S J063715-500414 and ACT-S J065207-551605, are not cross- 
identified with other catalogs, although we expect some of these 
to be real DSFGs, a hypothesis which will be tested through future 
follow-up observations. 

In order to estimate completeness of the catalog for each band, 
we added 100 synthetic sources of a single flux density to the 
455 deg 2 region of the ACT maps used in this study. The map filter- 
ing and source detection algorithm (Section 3) were run, and the re- 
sulting catalog checked for inclusion of the input sources. This pro- 
cedure was repeated in each band in intervals of 10 mJy in synthetic 
source flux density in the range 10-100 mJy. The completeness at 
each flux density was interpolated to estimate the completeness in 
flux density bins given in Table 2. 


5 COMPARISON TO OTHER CATALOGS 

To further characterize ACT sources, we consider cross- 
identifications with other catalogs. Given the number of ACT- 
detected sources and the ACT beam size, around 1 out of 10,000 
sources randomly placed in the 455 deg 2 area considered here 
would coincide with an ACT source. Figure 5 shows the ACT 
sources located in flux density space, and cross-identifications from 
some overlapping catalogs with extragalactic source detections in 
the southern hemisphere discussed in this section. The synchrotron- 
dominated sources occupy points near the power law S„ oc v a 
with spectral index a = -0.6. Dust-dominated sources will follow 
a line closer to a = 3.7, and are generally associated with sources 
detected with S/N > 5 at 218 GHz only. 

As our cross identifications come from radio catalogs at low 
frequencies, and dust-dominated sources lack strong synchrotron 
radiation, we expect fewer of the dust-dominated sources to be 
co-identified with a previously known radio source. Overlapping 
large-area, sensitive far IR catalogs would be the equivalent for our 
dust-dominated sources, but these do not exist. Even if they did, 
they would be too confused for positive cross-identification in this 
manner. At present, these unidentified sources are the subject of 



S i4s[ mj y] 

Figure 5. The ACT sources cross-identified in flux density space. Sources 
have been selected with S/N >5.0 in at least one frequency band (dotted 
lines). See text for a more thorough description of cross-identification statis- 
tics with other catalogs. The diagonal lines follow the power law S v oc u a 
with oc- 1.66 for the dashed line, distinguishing the o ' —0.6 (dash-dot) 
synchrotron-dominated source population from the a = 3.7 (dash-dot-dot) 
dust-dominated source population. 

individual follow-up studies (e.g., Ikarashi et al. 2011; Scott et al. 
2011; Lupu et al. 2012; WeiB et al. 2013). 

5.1 General Statistics of Identifications 

Of our 191 sources, many have cross-identifications with 
several catalogs. One hundred seventy-four are identified with 
sources in the 0.84 GHz Sydney University Molonglo Sky Sur- 
vey (SUMSS; Mauch et al. 2003), and 122 of those cross-identified 
sources also belong to the 4.85 GHz Parkes-MIT-NRAO (PMN) 
survey radio catalog (Wright etal. 1994). Fourteen sources were 
cross-identified with the Infrared Astronomical Satellite (IRAS; 
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Helou et al. 1988) at 12-100 fim. Sources were also identified with 
one or both of the Australian Telescope 20 GHz survey (AT20G; 
Murphy et al. 2010) and the 1. 4/2.0 mm SPT (Vieira etal. 2010) 
catalogs. Two of the remaining 8 unmatched sources were observed 
during the November 2010 ATCA follow-up campaign (Sec- 
tion 2.3). This leaves 6 sources (ACT-S J023600-530237, ACT- 
S J024430-541605, ACT-S J035034-524801, ACT-S J062747- 
512614, ACT-S J063715-500414 and ACT-S J065207-551605) 
with no identification in any of the available catalogs, all but one of 
which have dust-dominated spectra. The single unidentified source 
with a synchrotron spectrum, ACT-S J023 600-5 3023 7, falls near 
our S/N = 5 threshold at 148 GHz and could be a false detection. 

A NASA/IPAC Extragalactic Database (NED) search using a 
radius of 2' cross-identified 3 sources from the ACT catalog with 
galaxy clusters: Abell S0250, Abell 3391 and Sersic 037/01. In- 
creasing the search radius to 3' adds Abell 3128 and Abell 3395 to 
this list. Simulations of the microwave sky suggest that only « 3% 
of galaxy clusters have their 148 GHz signal contaminated at 20% 
or more (Sehgal et al. 2010). 


5.2 Matches to AT20G 

The AT20G survey, carried out from 2004 to 2008 at 20 GHz 
with follow-up at 5 and 8 GHz, covered 6.1 sr of the Southern sky 
to a flux density limit of 40mJy at 20 GHz (Murphy et al. 2010). 
Of our sources, 115 match sources in the AT20G catalog, listed in 
Table 4, all of which we classify as synchrotron-dominated. Given 
the AT20G completeness limit and the predominance of flat spectra 
for our 148 GHz-selected sources (Section 7), faint ACT sources 
may not have matches in AT20G. We received time on the Australia 
Telescope Compact Array to measure flux densities for sources in 
the earlier Ml 1 catalog that did not appear in AT20G, the results of 
which are described in Sections 2.3 and 7. 


5.3 Matches to the Planck Compact Source Catalogs 

The Planck satellite team has released two sets of all-sky cata- 
logs in nine frequency channels, including two very close to those 
used by ACT. The first set of catalogs made up the Early Release 
Compact Source Catalog (ERCSC; Planck Collaboration XIII 
2011). Very recently, a deeper catalog based on 1.6 complete 
Planck surveys has been released, the Planck Catalog of Compact 
Sources (PCCS; Planck Collaboration XXVIII 2013). At 143 and 
217 GHz, the Planck beam sizes are approximately 1' and 5' re- 
spectively. Like WMAP, Planck makes great gains by observing 
the entire sky simultaneously at several frequencies, thus captur- 
ing all bright sources. On the other hand, due to their large beam 
sizes, these experiments detect only the brightest sources. In the 
case of the PCCS, the detection threshold at 143 and 217 GHz is 
« 400 mJy; at lower flux densities the catalog completeness drops 
well below 80%. In the area treated in this study, Planck detects 3 
of the brightest ACT sources. 

For a more complete comparison with the PCCS, we exam- 
ined all bright ACT sources, combining the Equatorial and South- 
ern regions mapped by ACT. Almost 50 matches were found in 
this wider comparison. At 148 GHz, the agreement of ACT flux 
densities and the PCCS flux densities, properly color-corrected and 
(slightly) extrapolated to match ACT’s central frequency, is excel- 
lent, at the 1-2% level, despite scatter introduced by source variabil- 
ity. At 218 GHz, source variability plays at least as large a role, and 
PCCS flux densities are fa 5% higher than ACT’s. If we remove a 


couple of the most variable sources, however, the agreement of the 
218 GHz flux density scales improves to fa 1%. At both frequen- 
cies, Planck fluxes are on average slightly higher. Agreement at the 
fa 1-2% level, however, is well within the expected uncertainty in 
the ACT flux densities (Section 2). This agreement suggests that 
our flux density values are free of systematic error. 

5.4 Matches to the SPT Catalog 

The SPT study of compact sources (Vieira et al. 2010) was 
based on observations of a square patch of sky of 87 deg 2 centered 
at 05 h right ascension, having only fractional ACT overlap. Nev- 
ertheless, 2,304 of the 3,496 SPT candidate sources (those with 
S/N > 3 and flux densities > 4.4 mJy) fall within the ACT survey 
region. We find 32 cross-identifications with ACT sources when 
we search the SPT 1.4 mm (220 GHz) and 2.0 mm (145 GHz) cata- 
logs. Twenty-five of these were 148 GHz-selected and categorized 
in Vieira etal. (2010) as synchrotron-dominated. The remaining 
seven were detected by ACT with S/N > 5 only at 218 GHz; three 
of the seven also have SUMSS or IRAS cross-identifications. To 
make a direct comparison of the flux densities of ACT-SPT cross- 
identified sources, we used our derived flux densities, but without 
introducing two angular degrees of freedom in the deboosting (as 
in Vanderlinde et al. 2010), as these were not incorporated by the 
SPT analysis (Crawford et al. 2010). We find an overall offset at 
the 7 ±5% level at 148 GHz and 3 ±7% at 218 GHz, with ACT 
flux densities typically higher than those of SPT. At low flux, the 
ACT and SPT measurements of individual sources agree within 
their errors. However, at fluxes above 30mJy, individual sources 
may disagree at several sigma. We understand the discrepancy to 
be due to variability in radio sources, which is deferred to a future 
publication. 

Recently, using spectroscopic follow-up observations with 
the Atacama Large Millimeter/ submillimeter Array (ALMA), 
Vieira etal. (2013) and WeiB et al. (2013), derived robust red- 
shifts for about 18 of the SPT-detected dust-dominated sources 
from a larger sky area than in Vieira etal. (2010). One of these 
sources, found to be at redshift 5.66, matches ACT-S J034640- 
520505 which otherwise had no cross-identification. A second pre- 
viously unmatched source, ACT-S J002707-500713, was imaged 
with ALMA at 870 /^m, is resolved at the 0.5" scale, likely due 
to gravitational lensing, but awaits spectroscopic follow-up. Four 
other ACT dust-dominated sources that were cross-identified with 
the Vieira et al. (2010) catalog now also have robust redshifts. ACT- 
S J055 139-505800 is at redshift 2.123, ACT-S J053250-504709 is 
at redshift 3.399, ACT-S J052903-543650 is a source at redshift 
3.369, with a lens at redshift 0.140, and ACT-S J053817-503058 
has a redshift of 2.782, and a lens at a redshift of 0.404. It is likely 
that the few ACT dust-dominated sources which lie outside the SPT 
footprint are also lensed DSFGs. 


5.5 Matches to the IRAS Catalog 

Given the ACT beam size, a normal galaxy will be unresolved 
at redshifts z /S 0.05 or distances greater than 200 Mpc. Conse- 
quently, only very nearby objects appear extended in our maps. 
The IRAS source population consists primarily of local (z <C 1) 
dust-dominated ULIRGs; we thus expect our resolved sources to 
coincide with IRAS sources. We find 14 sources cross-identified 
with IRAS sources, out of a possible 829 IRAS sources that 
lie within the ACT survey area treated here. Two of these have 
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synchrotron-dominated spectra, both of which are nearby galax- 
ies. The twelve that show spectra dominated by dust re-emission 
are Galactic sources, in the Magellanic clouds, or are known 
nearby star-forming galaxies. For example, for two local resolved 
galaxies particularly bright at ACT frequencies, NGC1566 (ACT- 
S J04 1959-545622), a Seyfert two-arm spiral, and IC1954 (ACT-S 
J033 133-5 15352), a late-type spiral with a short central bar, ACT 
observes a higher 218 GHz flux density than at 148 GHz or 20 GHz, 
confirming the dust-dominated nature of their spectra. 

5.6 Dust-Dominated ACT Sources 

Of the 24 ACT dust-dominated sources in the catalog pre- 
sented here, 18 of these are cross- identified with either IRAS 
(12 sources) and/or SPT (9 sources, including the recent ALMA 
observations). ACT-S J05 1506-534420, ACT-S J05331 1-523827, 
and ACT-S J055 11 5-533435 were cross-identified with both cata- 
logs. Ten of the dust-dominated sources are cross-identified with 
sources in the SUMSS catalog, bringing the total for cross- 
identified sources to 19. The remaining 5 sources, ACT-S J024430- 
541605, ACT-S J035034-524801, ACT-S J062747-512614, ACT- 
S J063715-500414 and ACT-S J065207-551605, have no match- 
ing counterpart, and signal-to-noise at 218 GHz in the range 
5.27 < S/N < 6.35. Purity tests suggest that a couple of these could 
be spurious. However, given that three S/N « 6 sources are cross- 
identified with SPT, some of these sources are likely real detec- 
tions, high redshift galaxies lensed by intervening structure. Future 
follow-up observations are planned to clarify the nature of these 
detections. 

The main and supplementary 500 /j,m samples from the 
Herschel HerMES survey have lensing candidate densities 
of 0.14±0.04deg~ 2 and 0.31 ± 0.06 deg -2 , respectively 
(Wardlow et al. 2013). An analogous sample from Herschel ’ s 
H-ATLAS survey is 0.35 ± 0.16 deg -2 (Negrello et al. 2010). In 
the same bands as presented here, SPT finds 0.25 ±0.02 deg~ 2 
lensing candidates (Mocanu et al. 2013). This suggests that « 100 
such sources exist in the field here, of which we see the high flux 
tail. 


6 SOURCE COUNTS 

The completeness-corrected differential number counts for 
ACT sources based on the data in Table 2 are plotted in Figure 
6. For comparison, completeness-corrected number counts from 
Vieira et al. (2010) and Planck Collaboration XIII (201 1) are plot- 
ted as well. The ACT data fill in the flux density gap between the 
SPT and Planck catalogs at these frequencies, caused by the differ- 
ences between experiments in sky coverage and sensitivity to point 
sources. The effect of calibration error is to shift the flux bins by 
± 2% (± 2.4%) at 148 GHz (218 GHz; Section 2). 

The combined ACT, SPT and Planck total counts at 148 GHz 
and synchrotron-dominated source counts at 218 GHz are best fit 
by the recently developed Tuccietal. (2011) C2Ex radio source 
model. This model divides the ~ flat- or steep-spectrum blazar pop- 
ulation into BL Lac objects, with a region close to the AGN core 
dominating the observed emission (^ 0.3 pc), and flat spectrum ra- 
dio quasars (FSRQs) with a break frequency indicating that the 
source of emission arises from an emitting region further from the 
AGN core (0.3 to 10 pc). The de Zotti et al. (2005) model for radio 
source counts, while consistent with the counts below 0.1 Jy, over- 
predicts the counts at higher flux densities. Planck analysis finds 


that the de Zotti et al. (2005) model is consistent with their counts 
at frequencies up to 100 GHz, but over-predicts the counts at higher 
frequencies in the flux density region of « 1 Jy, though they begin 
to suffer from incompleteness below « 1 Jy. Our few « 1 Jy bright- 
est sources appear to be consistent with this finding, but lend little 
statistical significance. 

As well as the models for counts of radio sources, the source 
counts predicted for dusty starburst galaxies from Toffolatti et al. 
(1998), Limaetal. (2010), Bethermin et al. (2011) and Cai et al. 
(2013) are shown in Figure 6. The brightest infrared sources in 
the Toffolatti et al. (1998) model are 10 mJy. Given that all sources 
in the ACT catalog have flux densities greater than lOmJy, these 
models predict that the 148 GHz- and 218 GHz-selected ACT cata- 
logs should have few or no detected infrared sources. However, we 
find 24 sources that have dust-dominated spectra with flux densi- 
ties above 20 mJy, and Vieira et al. (2010) find 36 dust-dominated 
sources with flux densities above 15mJy. While Toffolatti et al. 
(1998) does not incorporate the effects of lensing on observed 
number counts, the models of Lima et al. (2010), Bethermin et al. 
(2011), and Cai et al. (2013) modelled dust-dominated galaxies to- 
gether with lensed high-redshift dust-dominated galaxies. How- 
ever, the small number of dust-dominated sources present in Figure 
6 lend little constraining power, and we leave it to future larger- 
area ACT studies to analyze the robustness of these models. The 
Bethermin et al. (201 1) model, falling roughly in the middle of this 
group of models, was used for our flux deboosting. 


7 SOURCE SPECTRA 

Source spectral energy distributions (SEDs) can be used to dif- 
ferentiate source types by their dominant emission mechanisms. 
Assuming the commonly used simple power law model S(u) oc 
v a , a negative a is indicative of sources dominated by synchrotron 
emission, such as radio galaxies. Sources with free-free emission 
dominating will have an index close to 0. The high-redshift SMG 
population will have spectra dominated by re-emission of their 
prodigious optical and UV flux by the surrounding dust in a grey- 
body spectrum, with indices expected to be greater than 2 and more 
typically 3-4. 

We can divide the ACT source population according to sev- 
eral broad spectral groups: classical steep (and steepening) spec- 
trum sources, sources that peak within the frequency range under 
consideration, and sources that show flat, rising or upturned spec- 
tra. 

The ACT catalog is dominated by synchrotron-dominated 
blazers, which have variable flux densities. This variation in flux 
density is not periodic. For any single source, then, inferences about 
its spectrum will depend on the epoch of observation, although not 
biased one way or the other. For the catalog as an ensemble, how- 
ever, a spectral study may give rise to insights about the average 
spectral behavior of the galaxy populations. 

Table 3 summarizes the median spectral indices between pairs 
of frequencies for various subsets of the ACT data. In obtaining 
the spectral indices we compare the deboosted flux densities from 
ACT with the raw flux densities from AT20G. The whole sample 
includes data from the 20 GHz ATCA follow-up observations taken 
in November 2010 (Section 2.3). 

Figure 7 shows the 5-20 GHz, 20-148 GHz and 148-218 GHz 
color-magnitude diagrams. The radio-selected AGN (black points) 
are predominantly characterized by steepening of the spectra, with 
only a few characterized by extremely flat or inverted spectra. There 
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Figure 6. Differential number counts of ACT-selected sources. Derived from Table 2 and corrected for completeness, the ACT differential source counts are 
plotted together with models of radio and infrared source populations. The Planck Collaboration VII (2013) data points and the SPT data points of Vieira et al. 
(2010) are also plotted, both of which are consistent with ACT counts. The ACT and Planck points have Poissonian errors (ler), whereas the SPT points 
include measurement and independent calibration errors. The flR models on this scale predicts number counts too low to be seen at 148 GHz. The ACT data 
are consistent with being dominated by radio sources at both frequencies. 


Table 3. Median Spectral Indices. 


Spectral Index 

Synchrotron (all) a 

Synchrotron (S 148 > 50 mJy) 

Synchrotron (Si 4 g < 50mjy) 

Dust-dominated 

£*5—20 

-0 15" 1 ” 0 ' 37 
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-0 07”*”°' 36 
uu/ -0.25 

-0 21“*” 0,33 
U - Z1 -0.41 
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-0 42“*" 0,32 
-0.26 
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u " 3 °— 0.29 

-0 43+ 0 - 34 
w ^ -0.33 


0148-218 

-0 55“*" 0 ' 60 
0.60 

-0 60"*”°‘ 20 
U.OU-0 20 

-0 5i+ 0 - 38 
U - J1 -0.71 

t 7+0.62 

^■'-o.se 


“Quoted errors are the 68% confidence levels of the distribution. 


is a clear trend towards more negative median spectral index with 
increasing frequency. Magenta points, showing rising spectra be- 
tween 148 and 218 GHz, denote sources ACT-classified as spec- 
trally dust dominated. They do not show up in the top two plots, 
indicating that these sources have flux densities falling below the 
detection threshold of the AT20G catalog. 

The average spectral indices between AT20G and ACT fre- 
quencies indicates an underlying source population made up of 
Flat-Spectrum Radio Quasars (FSRQ), a type of blazar, with AGN 
jet pointed along our line of sight (de Zotti et al. 2010). Ejected 
material flows through several shocked regions in the jet which lo- 
cally enhance the radiation (Marscher & Gear 1985; Valtaoja et al. 
1992). The observed spectral flatness is the superposition of 
many components with different turnover frequencies. At frequen- 
cies greater than approximately 100 GHz, however, Marriage et al. 
(2011), Vieira et al. (2010), and Planck Collaboration XIII (201 1) 
observe a steepening of the spectrum (ai 48 - 2 is ~-0.6) that un- 
til now had not been conclusively shown (Tucci et al. 2011). This 
change is possibly due to electron energy losses in the jet (“elec- 


tron ageing”) or the transition to the optically thin regime in the 
extended radio lobes. The underlying physical mechanisms have 
been contested for more than a decade and remain the subject of 
ongoing study, (e.g., Nieppola et al. 2008; Ghisellini & Tavecchio 
2008; Sambruna et al. 2010). 

All but one of the synchrotron-dominated ACT detections 
with flux density > 50 mJy have cross-identifications in AT20G. 
Below this flux density, the mean spectral indices of the population 
of ACT-AT20G cross-identified sources are biased toward the neg- 
ative by the incompleteness of the AT20G catalog below 100 mJy 
at 5 GHz and below 40 mJy at 20 GHz. It is therefore illustrative to 
further divide this subpopulation according to 148 GHz flux den- 
sity. 

Figure 8, a radio color-magnitude diagram, plots spectral in- 
dex 020-148 against 148 GHz source flux density. Black and grey 
points identify the Mil sources that were cross-identified with the 
AT20G catalog. Black points are for sources with flux densities 
> 50mJy, which represented a complete sample. Grey points de- 
note the fainter, incomplete sample. The blue points were obtained 
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Figure 7. Color-magnitude diagrams comparing the 5-20 GHz (top), 20- 
148 GHz (middle), and 148-218 GHz (bottom) spectral indices for ACT- 
ATCA cross-identified sources. The synchrotron-dominated radio galaxy 
population is dominated by sources which have consistently falling SEDs 
towards higher frequencies. 


by calculating the spectral index for the previously unmatched, 
lower flux density suhsample, followed up with ATCA (see Ta- 
ble 4 for AT2G0 source IDs with an asterisk). The population rep- 
resented by the blue points fills in the picture remarkably for the 
fainter (below 50m Jy flux density at 148 GHz, below 40mJy at 
20 GHz) population, especially in the region with spectral index 
of approximately zero. 

The unbiased Si 4 s >50mJy sample has a 20-148 GHz me- 
dian spectral index of -0.36j/Q2g. For the 5i48 < 50 mJy subsample 
prior to follow-up with ATCA (grey points only), the 20-148 GHz 
spectral index was (*20-148 = -0.52+q; 28- However, with the addi- 
tion of the synchrotron-classified ACT-selected sources with 2010 
ATCA follow-up data (grey and blue points), the full sample has 
20-148 GHz spectral index (*20-148 = - 0 . 42 to; 2 g. This supports 
the hypothesis that the lower flux contingent is probing the same 
population of synchrotron-dominated sources (blazers), and that 
there isn’t much evolution of their spectral index with flux. The 
full dataset suggests that the increased scatter at lower fluxes is due 
to variability (which will have a larger relative effect on the fluxes) 
and decreased S/N from flux error. 

For our dust-dominated sources, we find a median spec- 
tral index of (*i48-2is =3.72o'.86- F° r sources detected above 
5a at both 150 and 220 GHz, SPT derives 0150-220 =3.35 ±0.7 
(Mocanu et al. 2013). These indices are consistent with the find- 
ings of Dunkley et al. (2011), where the effective index of unre- 
solved IR emission as determined from ACT data is 3.69 ±0.14, 
and the best-fit mean spectral indices of afso -220 =3.86 ±0.23 
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Figure 8. Radio color-magnitude diagram using 20-148 GHz spectral in- 
dices for ACT-AT20G cross-identified sources. The data are divided be- 
tween flux densities at 50mJy at 148 GHz. The low flux sample was in- 
complete and suffered from selection bias that favored sources with more 
negative spectral indices. Data from the ATCA follow-up study for the Mil 
ACT 148 GHz sources not cross-identified with AT20G are shown as green 
squares. The high flux sample, denoted with black points, has a median 
spectral index -0.37. Prior to ATCA follow-up observations, the lower flux 
ACT sources cross-identified with the AT20G catalog (grey crosses) had 
a median spectral index of -0.52. Including the follow-up data, the more 
complete lower flux sample has a median spectral index of -0.43. 
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for Poisson-distributed sources, and (* 150-220 = 3.80 ±1.3 for the 
clustered component, found below the SPT detection threshold 
(Hall et al. 2010), a population expected to be dominated by dust- 
dominated sources. If we model dust emission as a modified black- 
body, in the Rayleigh-Jeans (RJ) limit the dust emissivity index is 
related to the effective spectral index as /3 = a — 2 = 1 .7, consistent 
with models of « 30K dust made from graphite and silicate grains 
(Draine & Lee 1984). We leave a more rigorous analysis involving 
redshifted greybodies, where the RJ limit is not as good an approx- 
imation, and the implications for star formation, to future work. 


8 CONCLUSION 

We have described the extragalactic source population at 148 
and 218 GHz found in a 455 deg 2 region of the ACT 2008 South- 
ern strip, centered on declination -52.5°. This updates the ACT 
148 GHz catalog by using the data released with Dunneretal. 
(2013), extends the results in Marriage etal. (2011) to two bands, 
and treats noise as more local, which in turn yields a higher S/N. 
We detect 191 sources above S/N = 5 in at least one of the ACT 
148 and 218 GHz frequency bands, spanning flux densities 14— 
1700 mJy (Table 4). Known redshifts are as high as « 6, with mea- 
surements ongoing. The catalog is estimated to be 100% pure and 
96.8% complete above 30mJy at 148 GHz, and 97.7% pure and 
80.2% complete above 30 mJy at 218 GHz. We have confirmed flux 
recovery of the pipeline, and jointly deboosted the flux densities in 
both bands. 

The multifrequency nature of our observations allows for in- 
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temal classification of sources into three broad classes of sources 
based on their spectra: synchrotron-dominated sources (the vast 
majority of which are cross-identified with radio catalogs), low- 
redshift dust-dominated sources with IRAS counterparts (typically 
ULIRGs), and dust-dominated sources also observed by SPT or 
with no cross-identification, either shown or assumed to be high 
redshift star-forming galaxies. This last class, only recently ob- 
served at millimeter wavelengths, has many of the properties ex- 
pected of the progenitor population of massive, modern-day, el- 
liptical galaxies, background SMGs whose flux has been magni- 
fied through gravitational lensing by a foreground galaxy or galaxy 
group. This interpretation, bolstered by population synthesis anal- 
yses (e.g. Thomas et al. 2005), is being validated with follow-up 
observations. 

A comparison with other catalogs shows that 97% of ACT- 
detected sources correspond to sources detected at lower or higher 
frequencies. The 148 GHz source counts are fit reasonably well 
by the C2Ex radio model of Tucci et al. (2011), the most cur- 
rent model for radio sources. According to the analysis of the 
average spectral indices derived from the combined AT20G and 
ACT datasets, the ACT data support the case for a spectral steep- 
ening toward higher frequencies above 100 GHz for AGN. The 
ACT dust-dominated source population has a median spectral in- 
dex, ai 48 - 2 i 8 , of3.7tg;!£. Properly linking these sources into the 
broader context of galaxy formation and evolution is of cosmolog- 
ical interest, and a goal of future work. 

The analysis presented here uses only the 2008 data of ACT’s 
Southern strip, representing only l/6 th of the data ACT obtained 
between 2007 and 2010. In future work, we will extend our anal- 
ysis of the source population to include the full dataset integrating 
both the southern and equatorial regions observed by ACT. The 
ACT Equatorial strip overlaps with deep Sloan Digital Sky Survey 
Stripe 82 observations (Annis et al. 2011). Thus as well as increas- 
ing the sky coverage and number counts for the ACT sources, future 
work (Gralla et al. in prep.) will include joint analyses with optical 
data. 
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